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ABSTRACT 

We investigate the IR spectral features of a sample of D-type symbiotic stars 
in order to constrain the emitting properties of coupled dust-gas particles across 
the whole system. In particular, analyzing unexploited ISO-SWS data, deriving 
the basic observational parameters of dust bands, and comparing them with re- 
spect to those observed in other astronomical sources, we try to highlight the 
effect of environment on grain chemistry and physics. We find strong amorphous 
silicate emission bands at ~10 /im and ~18 /xm in a large fraction of the sam- 
ple. The analysis of the ~10 /im band, along with a direct comparison with 
several astronomical sources, reveals that silicate dust in symbiotic stars shows 
features between the characteristic circumstellar environments and the interstel- 
lar medium. This indicates an increasing reprocessing of grains in relation to 
specific symbiotic behavior of the objects. A correlation between the central 
wavelength of the ~10 and ~18 fim dust bands is found. By the modeling of IR 
spectral lines we investigate also dust grains conditions within the shocked neb- 
ulae. Both the unusual depletion values and the high sputtering efficiency might 
be explained by the formation of SiO molecular, which are known to be a very 
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reliable shock tracer. We conclude that the signature of dust chemical distur- 
bance due to symbiotic activity should be looked for in the outer, circumbinary, 
expanding shells where the environmental conditions for grain processing might 
be achieved. Symbiotic stars are thus attractive targets for new mid-infrared and 
mm observations. 

Subject headings: binaries: symbiotic - infrared: stars - stars: circumstellar mat- 
ter - stars: individual: R Aqr, CH Cyg, V627 Cas, HM Sge, V1016 Cyg, RR Tel, 
V835 Cen, Hl-36, BI Cm. 

1. Introduction 

Symbiotic systems (SS) are currently understood as interacting binaries composed of a 
compact star, generally but not necessarily a white dwarf (WD), a cool giant star, which is 
at the origin of dust formation and ejection episodes, and different emitting gas and dust 
nebulae. In the past years, theoretical models (Girard & Willson 1987, Kenny & Taylor 
2005) as well as observations (Nussbaumer et al. 1995) have definitively shown that both 
the hot and cool stars lose mass through strong stellar winds which collide within and outside 
the system, hence creating a complex network of wakes and shock fronts which result in a 
complicated structure of gas and dust nebulae (Nussbaumer 2000). 

Two main phases characterize the symbiotic activity: the quiescent phases, during which 
the hot component releases its energy approximately at a constant rate and spectral distri- 
bution, and the active phases, characterized by a significant change of the hot component 
radiation. Near infrared (NIR) photometry revealed the continuing presence of a cool vari- 
able star in these systems (Wallerstein et al. 1984 and references therein) both before and 
after an active phase, hence indicating that the eruption episodes do not significantly affect 
the cool component. Contemporary UV and IR observations showed that while the UV flux 
decreases during the quiescent phases, the IR remains virtually constant, again suggesting 
that the IR source is not directly involved in symbiotic activity (Friedjung et al. 1984). 
Nevertheless recent works, both theoretical and observational, highlit that cool stars in SS 
might behave quite differently with respect to the isolated ones. For instance, it has been 
suggested that both the red giants and the Mira variables in SS have systematically higher 
mass-loss rates than typical galactic giants and Mirae (Mikolajewska 2000; Mikolajewska et 
al. 2002), and that the cool giants in SS are likely to rotate much faster than the isolated 
ones, leading to an enhanced magnetic activity as a main result, and probably to inhomoge- 
neous mass-loss and asymmetric dust-ejection episodes (Soker 2002). 
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In the '70s, on the basis of near- infrared (NIR) colours, SS were classified in S and D 
types (Webster & Allen 1975) according to whether the cool star (S-type) or dust (D-type) 
dominates the 1-4 fim spectral range. A peculiar aspect of the symbiotic phenomenon in 
D types pertains therefore quite naturally to the dusty environment. These systems show 
a broad IR excess which, since the first IR surveys, has been attributed to emission from 
circumstellar dust. It was found that the dust excesses have colour temperatures near 1000 K 
(Feast et al. 1983, Whitelock 1987); nevertheless it has been difficult to explain the observed 
IR spectrum by a single temperature, and nowadays theoretical models and observations 
too have been confirming that several dust temperatures should be combined in order to 
reproduce the NIR-MIR data (Anandarao et al. 1988, Schild et al. 2001, Angeloni et al. 
2006a, hereafter Paper I, Angeloni et al. 2006b, hereafter Paper II). 

In previous investigations of the SED from SS (Paper I) we have analyzed the different 
contributions to the continuum, with special emphasis on emission by dust in the IR spectral 
range. The results revealed the presence of two dust shells in many systems, pointing out 
that multiple dust shells are common in D-type systems: the inner and hotter (~ 1000 K) 
shells appear circumstellar, with a typical size of 10 14 cm, while the outer and cooler (T=300 
K) with sizes of ~ 10 15 cm are likely to be circumbinary. Very recent observations by Biller 
et al. (2006) confirm that in some systems, such as CH Cyg and HM Sge, even the inner 
dust shell surrounds both the stars, interacting in symbiosis with the shocked nebulae. 

Since the last '60s smooth, broad bands at 10 /zm and 18 /zm have been observed in the 
spectra of some late-type, variable stars and then in a wider variety of galactic environments. 
They were soon associated with silicate dust grains and were assigned to the Si-0 stretching 
and bending modes respectively. Furthermore, although the silicates with various composi- 
tion all display these spectral features due to Si-0 vibrational modes, there are differences in 
the spectral appearance that allow to constrain the exact composition of the silicate grains, 
as well. After the ISO "crystalline revolution" (Jager et al. 1998) the whole dust profiles 
have been interpreted in terms of mixed crystalline and amorphous structures (e.g. Bowey 
& Adamson 2002, Kemper et al. 2005, Tamanai et al. 2006). 

Various attempts have been carrying on to classify the different O-rich dust features 
on the basis of the spectral profile, by both space-based facilities, like IRAS and ISO (e.g. 
Little-Marenin & Little 1990; Sloan & Price 1995; Posch et al. 1999) and ground-based 
ones, like UKIRT (Monnier et al. 1998; Speck et al. 2000, hereafter S00). Some studies 
highlighted also a clear dependence on the environment conditions: for instance, the central 
wavelength of the peak A c , of the ~10 fim band is known to decrease from >9.8 fimm the 
circumstellar and diffuse-medium environments to 9.6 /iin in the Taurus Molecular Clouds 
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(TMC), while the FWHM increases from 2.5 to 3.3 /xm (Pegourie & Papoular 1985). 
These changes in the profiles seem represent the effect of mixing grains from different sources 
and of chemical and physical processing on the older TMC and Young Stellar Object (YSO) 
silicates (Bowey & Adamson 2002, Table 2) and there have been even some works which 
attempted, on the basis of a long series of laboratory experiments, to relate the ratio of the 
10/im and 18 /im band strength with the relative age of silicate (Nuth & Hecht 1990). 
An accurate analysis of observational infrared spectra, discriminating the role played by the 
various environment parameters, is not trivial since not only mineralogical properties but 
also size, shape and agglomeration of the grains can have an influence on the shape and 
intensity of the distinctive spectral features. Thus further laboratory studies on the silicate 
dust bands are essential for an appropriate interpretation of the observed spectra. 

In the present work we analyze unexploited ISO data in order to compare the silicate 
dust features in D-type systems with those of "normal", isolated late- type stars and other 
galactic environments. We refer to the classification scheme proposed by S00 who, thanks to 
an infrared survey of more than 140 stars, identified similarities and differences between the 
features of various types of stars. As we are interested in Mira variables, that correspond to 
the cool component in D-type SS, which represent ~ 20% of the SS in the Belczyhski et al. 
(2000) atlas, we will focus on the three dust feature groups proposed by S00 for AGB stars: 
1) broad AGB, where the feature extends from 8 to about 12.5 /im with little structure; 2) 
broad+sil AGB, which consists of a broad feature with an emerging 9.7 /im silicate bump; 
and 3) silicate AGB, which is the "classic" 9.7 fim silicate feature. The silicate AGB group 
is split further into four subgroups, A, B, C, D, which show a sequence of slightly different 
shapes from broader to narrow silicate features (S00, Fig.l). 

Our aim is to study the environmental conditions in SS, both in the dust shells and in 
the shocked, ionized nebulae, in order to investigate whether the "symbiosis" might affect 
the dust behavior, leaving its marks on the spectral profile by changing the chemical and 
physical properties of the grains, as well as the shell morphology and location, as indicated 
by Paper I. The sample and data reduction are described in §2, while in §3 the results are 
presented and discussed. Concluding remarks follow in §4. 

2. The sample: observations and data reduction 

The sample has been selected by cross-checking Belczyhski et al. (2000) atlas with the 
ISO Archive, and then by including those D-types which showed the best quality observa- 
tions. It consists of the same 9 SS analyzed in Paper I: 8 D type, namely, R Aqr (RAq), 
V627 Cas (V6), HM Sge (HM), V1016 Cyg (VI), RR Tel (RR), V835 Cen (V8), Hl-36 
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(HI), BI Bru (BI), and 1 S-type, CH Cyg (CH), whose actual classification is however still 
controversial (see Paper I and references therein). 

We present the spectra taken with the Short Wavelength Spectrograph (hereafter SWS - 
de Graauw et al. 1996), used in full-grating scan mode (AOT 01) and covering the wavelength 
range between 2.38 and 45.2 /im. The spectra come from the "Uniform database of SWS 
2.4-45.4 micron spectra" within the Highly Processed Data Products (HPDP) section of the 
ISO Archive. The database presents a complete set of all valid SWS full-scan 2.4-45.4 /im 
spectra processed from the last stage of the pipeline software and renormalized as uniformly 
as possible, thus representing the most processed form available from the ISO archive (see 
Sloan et al. 2003 for details about the algorithm used to generate the database). However 
some instrumental artifacts are still present, and in specific ranges and cases the S/N level 
is not sufficient to derive reliable quantitative measurements: for instance the SWS band 3E 
at 27.5 - 29 /im is known for its mediocre performance and beyond that wavelengths the 
spectra quality gradually deteriorates (see Paper II). Anyway this does not affect the analisys 
of dust amorphous band at all, but it is reflected only in a inability to determine the very 
accurate values of the atomic line observed parameters. Considering the aim of the present 
work, we thus preferred to give only upper limits to the emission lines ratios, exploiting HM 
Sge as a test-case in the analisys of the emitting nebulae (§3.2.1). In the following of the 
present work we focus on the 7-27 /im spectral range. 

In the analysis of the dust bands at ~10 /im and ~18 /im we decided not to subtract 
any sort of "continuum", but to rely on the baseline defined by the task SPLOT within 
the NOAO-ONEDSPEC IRAF package software (see below) during the fitting procedure. 
Indeed we verified that, even when the continuum appeared moderately steep, the profile 
shapes did not change to the extent as to modify the basic band parameters (like A c and 
FWHM, the full width half maximum) and to prevent the classification of dust features. 
This has been confirmed a posteriori by the classification of the only star that the present 
sample share with the S00 one, i.e. R Aqr. In fact, in agreement whit S00, we classified it 
in the silicate AGB C group (see §3.1). 

Furthermore, the identification of the continuum to be subtracted is not trivial: S00, for 
instance, subtracted a 8 /im normalized 3000K blackbody representing a stellar photosphere, 
by explicitly assuming that in the 8-13 /im range the continuum is dominated by the stellar 
contribution. In fact, the modeling of ISO-SWS spectra (Paper I) has shown that in the 
abovementioned range the continuum emission is predominantly caused by reradiation of 
hot (~1000K) dust from the inner shells (r~10 14 cm). In the longer wavelength domain, 
i.e. around ~ 20 /im, the continuum becomes even more composite, being the sum of 
contributions from at least two dust shells at different temperatures. Any kind of continuum, 
without a model taking into account the real physical scenario, appears hence artificial and 
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almost meaningless. As our profiles did not change their shape substantially, we chose a 
direct approach to the dust band analysis. 

To measure the dust profile parameters we used two different approaches. The dust 
profiles are clearly asymmetric therefore, by definition, the central wavelength of the band 
will not coincide with the peak position: a gaussian fit is then not able in reproducing the 
exact, whole feature profile. Nevertheless, as we focused in the central wavelength of the 
peak we performed, by means of the NOAO-ONEDSPEC task within the IRAF package 
software, a gaussian fit of merely the upper portion of the band, fairly approximable by 
gaussian curve. Then, to check this kind of approach, we tried to follow the same method 
of Bowey & Adamson (2002), i.e. to estimate the central wavelength of the peak and the 
FWHM of the bands "by eye". A good agreement between the two methods was found, 
even if this translates in moderately larger error bars. Nonetheless it is very unlikely that 
some not accounted for physical effects could move the real central wavelength of the band 
peak as much as to shift it out of the error bars. It is clear that these uncertainties do not 
invalidate the results presented in the next section. As regards the atomic emission lines, 
we performed a gaussian fit by means of the above mentioned NOAO-ONEDSPEC task in 
order to measure their main observational parameters, namely the central wavelength and 
the FWHM. 



3. Results and discussion 

In Fig. CD we present the 7-22 /zm spectra for the SS of our sample. The symbiotic 
character of such systems can be inferred by noticing that on top of the composite stellar- 
dust continuum, and even overlapped on the dust broad bands, many intense forbidden 
emission lines are visible. 

In §3.1 we analyze and compare the amorphous silicate spectral profiles which dominate 
the wavelength range, and in §3.2, by exploiting HM Sge as a test-case, we model the IR 
line spectra in order to constrain the emitting properties of coupled dust-gas particles. In 
particular, we would like to investigate whether dust grains can survive within the shocked 
nebulae and how their chemical properties might be affected by the symbiotic environment 
whole. 
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3.1. Silicate spectral features 

The main characteristics of the spectra presented in Fig. [1] are the two broad, easily 
recognizable dust bands at 10 fim and 18 /im, arising from the vibrational modes of Si-0 
bound. These bands dominate the whole spectral range in 7 out of 9 SS: Hl-36, showing a 
flatter continuum whose shape is due to a cooler (700-800K instead of 1000K) dust (Allen 
1983; Paper I) and BI, where no traces of amorphous silicate bands are recognizable and 
the spectrum appears almost dust featureless; nonetheless, it is worth noticing that for this 
last object the poor S/N of the ISO spectrum doesn't allow to conclusively state about the 
actual presence of some dust features. 

3.1.1. The ~10 fj,m silicate profile 

The SS displaying the 10 fim band might be related, on the basis of the profile shape, to 
the S00 classification scheme for AGB stars. A direct comparison of our Fig. [2] with Fig. 1 of 
S00 shows that R Aqr and CH Cyg have a profile typical of the AGB silicate C group, HM 
Sge, V1016 Cyg and RR Tel belong to the AGB silicate D class, and V627 Cas and V835 Cen 
lie in the AGB broad feature group. No evidence for the ~13 fim feature sometime present 
in AGB stars is found. In Tab. [Hand Fig. [3] we investigate the effect of SS environments on 
the band profiles, plotting the FWHM vs. the central wavelength of the peak A c . SS appear 
in a specific locus of the diagram between the proper circumstellar objects and interstellar 
environments, being this a clue to an increasing reprocessing of grains across the systems 
(see below). 

It is worth noticing that the only two stars that exhibit the silicate C group profile (R 
Aqr and Ch Cyg) and that lie, in Fig. [3j nearest to the proper circumstellar environments 
(open circles) display just one dust shell (T^ ~ 1000K) and share similar characteristics 
both in the continuum SED (Paper I) and in the physical conditions of the ionized nebulae 
(§3.2.2). This might be just a coincidence: therefore we would like to further investigate 
the potential implications of such a link in the symbiotic system in a forthcoming, dedicated 
paper (Contini 2006, in preparation). 

The three stars (HM Sge, V1016 Cyg, and RR Tel) known to have a very complex network 
of wakes and shock fronts and many emission lines on top of the continuum, as evidence of 
strong contributions from the nebulae, display a ~10 jum profile belonging to the same group, 
the AGB silicate D. Even if the FWHM of the 10 fim band appears practically identical 
to that of typical stellar objects, the central wavelength of the peak is clearly at longer 
wavelengths (A c > 10.1 jttm), probably suggesting irregularity in shapes of dust grains (Min 
et al. 2003). These profiles are well fitted by a combination of Mg-rich silicates and porous 
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amorphous alumina and are interpreted in the light of a dust condensation sequence that 
can still be related to the specific condition of the isolated AGB star environments (Tielens 
1990, Sloan & Price 1998, Speck et al. 2000). 

As previously mentioned, V627 Cas and V835 Cen belong to the AGB broad feature group. 
In this case not only the A c appears at longer wavelength with respect to the other classes, 
but also the FWHM is larger (FWHM ~ 2.9-3.1 /xm). One of the interpretation of the larger 
width of the AGB broad feature group profile suggests a ~2 /im-wide excess centred at 11.5 
/im which has been attributed to the shifted peak of larger micrometer-sized grains (Jaeger 
et al. 1994). The presence of these larger grains could be a clue to crystalline silicates (van 
Diedenhoven et al. 2004), which have been observed in the D'-type SS HD330036 (Paper 
II). New high resolution and high S/N observations would be welcome considering the other 
species of silicates disclosed by the ISO "revolution". As a matter of fact these data might 
reveal the sharp substructures due to crystalline silicates which appear superimposed on 
the amorphous silicate broad features. A detailed analysis of the crystalline band would 
then be helpful in determining the physical parameters of these environments, thanks to 
derived estimates of mass-loss rates (Kemper et al. 2001), dust temperature (Bowey et 
al. 2002) and dust geometrical distribution (Molster et al. 2002a,b,c). On the other side, 
new laboratory measurements are necessary for a further verification of this insight under 
varying environmental parameters, allowing for a better correlation between the profile shape 
variations and the astrophysical condition parameters. 

To summarize, from this spectral comparison it follows on that dust in D-type SS 
appears, at least in some cases, different from the dusty environment of "normal" isolated 
late-type evolved stars. In particular, the sample seems to indicate a kind of sequence 
that moves towards objects with multiple and complicated structures of nebulae and shells. 
We suggest that for these systems (the ones that differ more from typical circumstellar 
environment on the basis of Fig. [3]) the different profile features can be interpreted in the 
light of mixed properties of the grains, being the grains responsible of the bulk of the emission 
both in the inner "circumstellar shell" and in the outer, circumbinary one (Paper I). The 
symbiotic behavior seems thus to play a dominant role not only in the dynamical parameters 
such as location and distribution of dust by affecting both the shell radii and shapes, but 
also in the chemical properties of the grains, which experience physical conditions both of 
circumstellar and ISM environments. The signature of dust chemical disturbance due to 
symbiotic activity should then be assigned to the outer shells of the whole system, where the 
environmental conditions for grain processing and even crystallization can, in some cases, be 
achieved (Paper II). 
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3.1.2. The ~18 jum silicate profile 

Fig. [I] shows that several SS have a second bump at ~ 18 /im, usually attributed to 
Si-0 bending mode in amorphous silicates. Even in this case, many emission lines are seen 
on top of the continuum and of the dust band. 

The strength of the ~18 /im profile decreases along with that of the ~10 /xm one; as these 
two features arise from different vibrational modes of the same chemical bound (Si-O), we 
wonder if the dust bands answer in a somewhat correlated ways to the same environmental 
conditions, probably periodically varying in reaction to the Mira periodic pulsations. 
To verify this idea we plot in Fig. S] (see also Tab. [T]) the 10 /im vs. 18 /im A c . Interestingly, 
we find a clear increasing correlation between these two observational parameters: the longer 
the A c (10), the longer the A c (18). The least square method was performed in order to find 
the best linear fit to the data: A c (18) = 1.39A C (10) + 3.73. 

Unfortunately the sample is poor in number of objects, and limited to D-type SS. We 
might expect that a similar behavior of amorphous silicate dust bands could also be present 
both in other circumstellar envelopes and galactic environments. It would be interesting 
if, after verified and constrained on a larger sample of different astrophysical environments, 
this relationship will gain some chemical and physical significances on the basis of further 
laboratory investigations. 

3.2. Dust and gas conditions in the ionized nebulae 

In SS two main shocks develop from collision of the winds from the stars, one propagat- 
ing outwards the system, and the other propagating in reverse towards the WD, between the 
stars. The nebulae downstream of the shock fronts are heated and ionized by the photoion- 
izing flux from the hot star and by shocks. In order to include this theoretical framework in 
our model, we complete our investigation on dust features in SS environments by calculat- 
ing dust and gas conditions downstream of the shocked nebulae. We will constrain the gas 
and dust physical conditions by modeling the strongest emission line ratios. The IR spectra 
show lines from a large range of ionization levels, including [FeVII] and [SiVII]. Radiation 
from the WD cannot lead to gas temperatures > 2-3 10 4 K in the nebulae, while the colli- 
sion processes heat the gas to temperatures which depend on the shock velocity (Toe V|). 
The observed lines are forbidden, therefore hardly emitted from the shells where very high 
electron densities lead to collisional deexcitation. Moreover, the high ionization level lines 
indicate high shock velocities characteristic of the nebula downstream of the reverse shock. 

The models which best fit the line ratios are calculated by SUMA (Viegas & Contini 
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1994; Contini 1997), a code that simulates the physical conditions of an emitting gaseous 
cloud under the coupled effect of ionization from an external radiation source and shocks, 
and in which line and continuum emission from gas are calculated consistently with dust 
reprocessed radiation as well as with grain heating and sputtering processes. Its calculation 
models have been already applied to several symbiotic stars, e.g. AG Peg (Contini 1997, 
2003), HM Sge (Formiggini, Contini & Leibowitz 1995), RR Tel (Contini & Formiggini 
1999), He2-104 (Contini & Formiggini 2001), R Aqr (Contini & Formiggini 2003), as well as 
to nova stars (V1974, Contini et al. 1997 - T Pyx, Contini & Prialnik 1997) and supernova 
remnants (e.g. Kepler's SNR, Contini 2004). The composite models are characterized by 
the parameters referring to the radiation flux (the temperature of the hot star, T*, and 
the ionization parameter, U), by those representing the shock (the shock velocity, V s , the 
preshock density, n , and the magnetic field, B ), and by those related with dust (the dust- 
to-gas ratio, d/g and the initial grain radius, a gr ). The relative abundances to H of He, C, 
N, O, Ne, Mg, Si, S, Ar and Fe are also accounted for. 



3.2.1. HM Sge: a test case 



We refer to the emission line measurements of SWS spectra for HM Sge by Schild et al 
(2001). In Table [2] the observed line ratios to [NeVI] are compared with model results. We 
consider first the Ne lines, [NeVI]7.65, [NeV]14.32, [NeIII]15.55, and [NeII]12.82, representing 
four different ionization levels of the same element. [NeVI] is the strongest line while [Nell] is 
very weak, indicating that the temperature of the gas is high enough across the whole nebula 
to prevent recombination. This constrains the model to matter-bound, suggesting that the 
emitting nebula is located between the stars, downstream of the shock which propagates in 
reverse towards the hot component. The relatively high velocities indicated by the FWHM 
of the lines (~ 500 kms -1 , Table [2]) confirm that we are dealing with the reverse shock. 
Such high velocities are in agreement with a plasma temperature of 10 6 6 K deduced from 
the X-ray observations by Miirset et al. (1997). The profile of the density and temperature 
downstream which result from the hydrodynamical equations determine the stratification of 
the ions and thus the intensity ratios of the lines. The best fit to the data is obtained with 
a preshock density no=5.10 5 cm -3 which leads to densities downstream n ~ 10 8 cm -3 , and 
with T* = 160 000 K and U=l. The corresponding ionizing photon flux is F u = 6 10 26 cm -2 
s _1 . If the radius of the ionizing source is Hwd ~ 10 9 cm and the distance to the nebula is 
r, we have F u K^ VD /r 2 = U n c, where c is the velocity of light. We obtain r=7.8 10 13 cm. 
The geometrical thickness of the nebula, D= 1.4 10 14 cm, which is lower than the binary 
separation (3 10 14 cm, Richards et al. 1999) is found by modeling. 
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Once the Ne line ratios are reproduced by the model, the relative abundances of the 
elements result from the agreement of the calculated with the observed line ratios. However, 
the elements different from Ne are represented by only one line in the observed spectrum, so 
the relative abundances are less constrained. 

In Fig. [5] top, the profiles of the electron density, n e , of the electron temperature of the 
gas, T e , as well as of the temperatures of dust grains, T gr , corresponding to different initial 
grain radii are shown across the nebula. The vertical thin solid line shows that most of the 
recombination region is cut off. The fractional abundance of the Ne ions (top) and of the 
ions corresponding to the other observed lines (bottom) are shown in Fig. [6] as a function of 
the electron temperature. 

In Table [3] our results are compared with those of previous works. We find O/H and Ne/H 
higher than solar by factors of ~ 2.3 and 4.8, respectively, revealing the strong contribution 
of the WD wind to the nebula between the stars (Nussbaumer & Vogel 1990, Weidemann 
2003). Moreover this allows us to postulate that oxygen had an initial (i.e. pre-1975 outburst) 
abundance as high as to compensate any depletion caused by subsequent grain formation 
processes. The other metals such as Mg, Si and Fe, which are easily trapped into grains, 
are instead depleted relatively to the solar values. Interestingly, the depletion of Fe does not 
appear correlated with that of Mg and Si, being instead more than twice larger. If iron was 
incorporated in the silicate lattice it would return to the gas phase along with the magnesium 
and silicon when the silicate is destroyed (see below); iron then must be trapped in different 
dust material, maybe in the form of refractory species such as metallic iron or iron oxide, 
so that it can stay in a solid phase while the magnesium and silicon are returned to the gas 
phase. This insight, which provides a natural explanation of the unusual depletion pattern 
observed also in several lines of sight through the interstellar medium (Welty et al. 2001; 
Sofia et al. 2006) has been already suggested by Sofia et al. (2006) and consistently confirmed 
by a very recent theoretical study by Min et al. (2006). As a further and complementary 
proof we deduce our results on the basis of the analysis of atomic emission lines. 

The rate of erosion (sputtering) of the grains can be high for V s > 200 kms -1 , leading 
to a complete destruction of the grains. Notice that even the largest ones (a gr = 1 /im) are 
sputtered within a relative small distance from the shock front (Fig. top). Mg, Fe and 
particularly Si should therefore return to the gaseous phase before the temperature drop to 
the values suitable to the line emission. This indicates that other effects must be considered, 
e.g. the lifetime of grain survival which depends on the distance downstream covered by the 
grains (Shull 1978, Dwek 1981). Small grains are stopped after distances of ~ 7 10 12 cm, 
while large grains after ~ 3 10 13 cm (Contini et al. 2004) which are of about the same order 
of the sputtering distances shown in Fig. 

An intriguing issue related to grain processing is that the return of silicon to the gas 
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phase might take place in the molecular form (SiO, SiC^), hence explaining the apparent 
contradiction between the high efficiency of sputtering processes we deduce from the shocks 
and the observed depletion calculated by the atomic lines. 

It is known that SiO is a very reliable tracer of shocked gas in many astrophysical environ- 
ment, and its enhancement in the gas phase has been already observed both in the bipolar 
outflows of young stellar objects (Martin-Pintado et al. 1992; Schilke et al. 1997) and in 
galaxies (Garcia-Burillo et al. 2000, 2001b; Usero et al. 2006). Moreover, there is growing 
evidence that different moleculae trace distinctly different velocity regimes in shocks, with 
SiO associated with more energetic events, i.e., those potentially more efficient in processing 
dust grains (Garay et al. 2000). Although qualitatively similar, shocks characterized by dif- 
ferent velocity regimes are expected to develop to a different extent dust grains in molecular 
gas: fast shocks can destroy the grain cores, liberating refractory elements to the gas phase, 
while slow shocks could only be able to process the icy grain mantles. One can then conclude 
that an increase in the typical velocity regime of shocks will favor an enhancement of the 
abundance of SiO in the shocked gas. 

The shock velocity characteristic of these processes, i.e. the threshold value above which the 
dissociation of SiO takes place, are usually V s < 60-70 km s _1 (Mart in- Pint ado et al. 1992). 
At the bottom of Fig. [5] the profile of the shock velocity downstream which results from 
the solution of the continuity equations for our model is shown, along with the [SiVII] line 
formation region. As one can see, at distances sufficiently large from the shock front the 
velocity drops to V s < 100 kms" 1 , indicating suitable conditions for molecular shock chem- 
istry. 

The interpretative framework is still in progress under many aspect, and it needs observa- 
tional data to be more constrained. Unfortunately, so far, there is a complete lack of spectral 
observations for SS in the mm spectral region. We then call for mm observations at SiO 
transition frequencies, e.g. ~ 217 GHz J=5-4 (A=1.38 mm), ~ 130 GHz J=3-2 (A=2.3 mm), 
~ 86 GHz J=2-l (A=3.4 mm), in order to trace the shock spatial profile within and around 
the symbiotic nebulae and to reveal up to what degree the sputtering processes are at work. 

Consider now the nebula downstream of the shock front expanding outwards the system 
which is characterized by a higher D, lower velocities and lower densities, hence being the 
source of the rich optical and UV spectra reported and studied extensively by previous 
works (e.g. Nussbaumer & Walder 1993, Formiggini et al. 1995). The grains can survive to 
destructive processes and contribute to the IR emission even at larger wavelengths, because 
of the lower temperature (T gr ~ 300K, Paper I) of dust. 

The overall conditions of dust and grains within the outer expanding nebulae are approaching 
those of the ISM up to a complete merge, indicating that the signature of dust chemical 
disturbance due to symbiotic activity, i.e. the chemical and physical reprocessing of the 
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grains suggested on the basis of the spectral analysis, might take place and should be looked 
for in these circumbinary structures. 

3.2.2. The other objects 

Fig. [T] shows that some objects of the sample have very similar spectra, both in the 
dust and in the emission lines features, while some others lack highly ionized lines (e.g. R 
Aqr, CH Cyg) and/or dust band (e.g. BI Cru, Hl-36). 

For objects with an evident emission line spectrum the Ne spectral line ratios, which 
are independent from relative abundance effects, are compared in Table HI It can be noticed 
that the general trend is very similar to that found for HM Sge, except RR Tel which shows 
[NeIII]/[NeVI] lower by a factor < 2, and Hl-36, which displays a very bright [OIV]25.89/im 
emission line, as well as no evident dust spectral feature (see §3.1). In Table Owe present 
the FWHM measured for the brightest lines. The derived velocities indicate that also for 
these objects, as for HM Sge, we are dealing with the inner shocked nebula, which emits 
the detected atomic emission lines and where dust grains can hardly survive the harsh 
environmental conditions. Considerations similar to those reported in the previous section 
(§3.2.1) about the molecular emission are still applicable. 

It is noteworthy that R Aqr and CH Cyg, known to have just one hot (T gr ~ 1000K) 
and circumbinary dust shell, do not show emission lines from highly ionized elements (Paper 
I). A plausible explanation for the absence of such strong atomic features could lie in the 
high density of the inner shocked nebula, up to 10 9 cm -3 (Contini 2006, in preparation), 
that prevents the formation of forbidden emission lines. 

Unfortunately for some objects ISO spectra are too noisy to reveal fainter spectral 
feature (V835 Cen, BI Cru). Furthermore, there are still insufficient data in some spectral 
ranges for a consistent and thorough analysis, and even for concluding about the actual 
symbiotic nature of some stars, as V627 Cas (Belczyhski et al. 2000, Angeloni et al. 2006, 
in preparation). Thus we postpone to forthcoming papers a detailed discussion of the single 
objects. 

4. Concluding remarks 

We have investigated the IR spectral features of a sample of D-type symbiotic stars in 
order to constrain the emitting properties of coupled dust-gas particles across the whole sys- 
tems. Having extracted the basic observational parameters of both dust bands and atomic 
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emission lines from unexploited ISO data and having compared them with those from other 
astronomical sources, we have attempted to highlight the effect of the environment on grain 
chemistry and physics, in particular whether dust grains can survive within the shocked neb- 
ulae and how their chemical properties might be affected by several depletion and sputtering 
patterns. 

The main conclusions of this work are the following: 

• The strong amorphous silicate emission features at ~10 /im and ~18 /im are found 
in 7 out of 9 objects of the sample. A qualitative analysis of the ~10 /im dust band 
indicates that the profile shape is similar to that of "standard" isolated AGB stars. 
Nevertheless a quantitative comparison of A c vs. FWHM for the 10 /iin silicate profile 
with other astronomical sources shows that SS appear between the characteristic cir- 
cumstellar objects and interstellar environments, which is believed to be a clear clue 
of an increasing, ongoing reprocessing of grains. 

• A correlation between the central wavelength of the ~10 /xm and ~18 /im dust bands 
is found. It is likely that the dust bands answer in a somewhat correlated ways the 
same environmental conditions, probably periodically varying in response to the Mira 
periodic pulsations if the dust grains responsible of the bulk of the emission are mainly 
located in the hot, circumstellar shell. Both new high resolution and high S/N ob- 
servations, as well as new laboratory measurements under varying environmental pa- 
rameters, are necessary for a further verification of this insight, even thanks to larger 
samples. 

• The modeling of IR spectral lines allowed to constrain the emitting properties of cou- 
pled dust-gas particles within the shocked nebulae, showing that sputtering processes 
are efficiently at work in the inner nebulae where high shock velocities are reached. The 
signature of dust chemical disturbance due to symbiotic activity should then be looked 
for into the outer, expanding, circumbinary shells where the environmental conditions 
for grain processing and even crystallization might in some cases be achieved. 

• Both the unusual depletion values and the high sputtering efficiency could be explained 
by the overabundant presence of SiO moleculae, known to be a very reliable shock tracer, 
hence making SS interesting targets for mm observations. 
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Table 1: A c and FWHM of the 10 /im silicate profile, A c of the 18 /im silicate profile, and 
relative uncertainties [/im] for the objects in Figs. [3]and|4l 



Object 


A c (10) 


err A c (10) 


1 WliJVl(lO) 


err b WliM(10) 


A c (18j 


err A c (18) 


Symbiotic 














R Aqr 


10.0 


±0.1 


2.5 


±0.1 


17.5 


±0.2 


CH Cyg 


10.0 


±0.1 


2.5 


±0.1 


17.6 


±0.2 


V627 Cas 


10.4 


±0.1 


2.9 


±0.1 


18.2 


±0.2 


HM Sge 


10.3 


±0.1 


2.6 


±0.1 


17.9 


±0.2 


V1016 Cyg 


10.1 


±0.1 


2.7 


±0.1 


17.9 


±0.2 


RR Tel 


10.1 


±0.1 


2.6 


±0.1 


18.1 


±0.2 


V835 Cen 


10.4 


±0.1 


3.1 


±0.1 


18.3 


±0.2 


Novae 














Aql 82 


10.2 


±0.1 


2.9 


±0.1 






Her 91 


10.7 


±0.1 


3.0 


±0.1 






Cas 93 


9.6 


±0.1 


1.8 


±0.1 






Circumstellar 














/i Ceph 


9.9 


±0.05 


2.5 


±0.1 






CU Cep 


9.9 


±0.05 


2.5 


±0.1 






U Aur 


9.8 


±0.05 


2.8 


±0.1 






Diffuse medium 














Cyg OB2 n.12 


9.7 


±0.05 


2.6 


±0.1 






), Orion Trap., TMC 














HL Tau 


9.6 


±0.05 


3.2 


±0.1 






Orion Trap. 


9.6 


±0.05 


3.5 


±0.1 






Taurus-Elias 7 


9.6 


±0.05 


3.5 


±0.1 






Taurus-Elias 13 


9.5 


±0.05 


3.4 


±0.1 






Taurus-Elias 16 


9.6 


±0.05 


3.3 


±0.1 






Taurus-Elias 18 


9.5 


±0.05 


3.2 


±0.1 
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Table 2: emission lines, normalized to [NeVI] 7.65 jum, in the SWS spectrum of HM Sge. 



A [/im] 


ID 


iP [eV] 


FWHM [/H 


observed 


model 6 


4.53 


[ArVI] 


91.01 


0.01 


0.09 


0.09 


5.61 


[MgV] 


141.27 


0.04 


0.19 


0.18 


6.49 


[SiVII] 


246.52 


0.03 


<0.07 


0.07 


7.65 


[NeVI] c 


157.93 


0.03 


1 


1 


9.53 


[FeVII] 


124.98 


0.04 


0.05 


0.05 


10.51 


[SIV] 


47.22 


0.04 


0.04 


0.04 


12.81 


[Nell] 


40.96 


0.07 


<0.07 


0.001 


14.32 


[NeV] 


126.21 


0.06 


0.13 


0.12 


15.55 


[Nelll] 


63.45 


0.06 


0.08 


0.08 


24.31 


[NeV] 


126.21 


0.11 


0.06 


0.04 


25.89 


[OIV] 


77.41 


0.10 


0.09 


0.08 



a: adapted from Schild ct al. 2001; 

b: input parameters: V s =500km s _1 ; n =5 10 5 cm" 3 ; B =10" 3 G; d/g=4 10 ~ 5 (by mass); & gr =l, 0.5, 0.2 /urn; D=1.74 10 
cm; T gas =1.85 10 4 K; T*=1.6 10 5 K; U=l; 
c: absolute flux: 74 ± 0.6 10~ 12 erg cm -2 s _1 . 



Table 3: the comparison of relative abundances for HM Sge. 



element 


model 


NV90 


solar 


He/H 


0.11 


0.10 


0.085 


O/H 


1.5(-3) 


1.0(-3) 


6.6(-4) 


Ne/H 


4.0(-4) 


4.0(-4) 


8.3(-5) 


Mg/H 


1.5(-5) 


1.0(-4) 


2.6(-5) 


Si/H 


8.0(-6) 


2.0(-4) 


3.3(-5) 


S/H 


3.6(-5) 


1.0(-4) 


1.6(-5) 


Ar/H 


6.3(-6) 


6.0(-6) 


6.3(-6) 


Fe/H 


6.0(-6) 


3.0(-5) 


4.0(-5) 



NV90: Nussbaumer & Vogel (1990). 
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Table 4: the brightest Ne emission lines" visible on top of continuum. 



Symb 


[NeV] 14 . 32 


[NeIII] 15 .55 


[NeV] 2 4.3i 


R Aqr 








CH Cyg 








V627 Cas 








HM Sge fe 


0.13 


0.08 


0.06 


V1016 Cyg 


<0.3 


<0.1 


<0.2 


RR Tel 


<0.2 


<0.05 


<0.1 


BI Cru 








V835 Cen 


<0.05 






Hl-36 


<0.5 


<0.2 


<0.15 



a: normalized to [NeVI] emission line @ 7.65 ^m; 
b: adapted from Schild et al. 2001. 



Table 5: FWHM [/zm] of the brightest atomic emission lines. 



Symb 


[NeVI] 7 . 65 


[NeV] 14.32 


[Nelll] 15.55 


[NeV] 2 4.3i 


[OIV] 25 .89 


R Aqr 












CH Cyg 












V627 Cas 












HM Sge 


0.03 


0.06 


0.06 


0.11 


0.10 


V1016 Cyg 


0.03 


0.05 


0.05 


0.11 


0.11 


RR Tel 


0.03 


0.06 


0.06 


0.09 


0.10 


BI Cru 












V835 Cen 


0.03 


0.06 








Hl-36 


0.03 


0.05 


0.06 


0.10 


0.11 
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Fig. 1. — The 7-27 /zm region of the ISO-SWS spectra for the SS of the sample. Notice the 
similarities and the differences in the continuum SED, in the dust band profiles at ~ 10 and 
~ 18 (Km (e.g. likely absent in BI Cru, even if the poor S/N of the spectrum doesn't allow 
to conclusively state about it) and in the intensities of atomic emission lines. 
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Fig. 2. — The 7-15 /im range. The amorphous silicate profile at ~10 /jm is dominant in 7 
out of 9 SS. 
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Fig. 3. — Effect of environment on the 10 /im silicate profile. Filled circles: SS; open 
pentagon: Novae (Aql82, Her 91, Cas 93 - Evans et al. 1997); open circles: circumstellar 
(CU Cep, U Aur, \x Cephei; Russell et al. 1975); open star: diffuse medium (Cyg OB2 n. 
12); open triangles: YSO and TMC (HL Tau, Taurus-Elias 7 and 18) and Orion Trapezium 
(Gillett et al. 1975); open squares: quiescent TMC (Taurus-Elias 16 and 13). 
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Fig. 4. — The peak wavelength of the 10 /im vs. 18 /im amorphous silicate bands. The 
dashed lines is the best linear fit with the data. 
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Fig. 5. — Top: the profile of electron density, electron temperature and grain temperature 
relative to different initial grain size (in jum) downstream across the reverse shocked nebula 
in HM Sge. The vertical thin solid line represents the edge of the matter-bound model. 
Numbers near the grain temperature line indicate the initial grain radius, a flr . Bottom: the 
profile of shock velocity, V s , and the distribution of the Si +6 /Si fractional abundance (thick 
solid line). 




Fig. 6. — The fractional abundance of the ions downstream as function of the temperature 
across the nebula for HM Sge. 



